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Considerable insights can be obtained regarding the nature of highly excited states by computing 
the eigenbasis expectation values of the resonance operators associated with an effective spectro- 
scopic Hamiltonian. The expectation values are related to the parametric derivative of the eigenval- 
ues with respect to specific resonance strengths i.e., level velocities. Sensitivity of the level velocities 
to the existence of closed orbits in the underlying classical phase space provides for a dimensionality 
independent route to a dynamical assignment of the states. In this letter, taking ClfBrClF (polyad 
P = 5) and the [16, 0]®^ bend polyad of C2H2 as examples, we show that the level velocities can 
signal the birth of new modes and highlight sequences of localized eigenstates. 
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I. INTRODUCTION 

Dynamical assignment of the highly excited eigenstates 
of a polyatomic molecule is a topic of significant interest 
to the chemical physics community 2^, ^ M- Insights 
into the nature of the excited states provides for a bet- 
ter understanding of the phenomenon of intramolecular 
vibrational energy redistribution (IVR) occuring in the 
molecule m, |3j, y|. Assignment of the low energy vibra- 
tional states in terms of the usual normal mode quantum 
numbers is relatively straightforward. However with in- 
creasing energies the normal mode quantum numbers are 
no longer conserved due to perturbations that strongly 
mix the zeroth order modes and the molecular vibrational 
Hamiltonian'^ takes the form: 
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wherein Hq represents the Dunham expansion, Xij are the 
anharmonic constants and Vk are the various perturba- 
tions. The operators dk and represent the annhilation 
and creation operators for the fc**^ mode respectively. The 
normal mode quantum numbers Vi — a|ai, sufficient for 
assigning low energy eigenstates, do not commute with Vk 
and hence cannot be used to assign the eigenstates of H. 
The above Hamiltonian is usually generated from a fit to 
the experimental spectrum in the absence of a global ab 
initio potential energy surface. Even in the case when an 
accurate potential energy surface is at hand it is possible, 
and useful, to generate a Hamiltonian of the above form 
using the canonical Van-Vleck perturbation theory|^0|. 
The advantages of using Hamiltonians of the above form 
stems from the fact that the classical limit Hamiltonian 

if (I, 6) = ffo(I) + 2^ rfe/fe(I) cos(mfc • 6) (2) 
fc 

is easily obtained via the correspondence a <-> VTe^*^. 
The /-dimensional vector nik has integer components 
{ri,r2, . . . ,rf) and (I, 9) are the action-angle variables^ 



corresponding to Hq. The order of a resonance 1^ is de- 
fined as Ok = X]j=i I'^il- "^^^ perturbations Vk{l,0) = 
/fc(I)cos(mfc • 6) are called resonant because the con- 
dition nik ■ = implies a specific commensura- 
bility or locking between the unperturbed frequencies 
= dHo{T) /dl. Such resonances are responsible for 
energy flow through the molecule and lead to breakdown 
of the zeroth order, low energy quantum numbers. Clas- 
sical dynamics of resonant Hamiltonians is historically 
very rich and have been studied in great detail0. In 
particular a detailed understanding of the classical dy- 
namics of H{I, 6) is essential for any dynamical assign- 
ment of the quantum eigenstates of H. Furthermore the 
effective Hamiltonian approach utilizes the state space 
perspective which offers considerable advantages towards 
understanding IVRj^]. 

Clearly any assignment of the eigenstates of H for- 
mally requires the existence of a sufficiently large set of 
good quantum numbers. Such a set does not exist in 
general and hence assignment implies the existence of at 
least approximate or quasi quantum numbers. By neces- 
sity such approximate quantum numbers are conserved 
for a certain time period dictated by the dynamics of the 
system. Support for the notion of quasi quantum num- 
bers partially comes from the experimental observation^^ 
of hierarchical IVR in molecules and the fact that most 
molecules are not ergodic even at fairly high energies 
In this sense deciphering the quantum numbers, exact 
or approximate, invariably implies the knowledge of the 
underlying dynamics. Consequently one speaks of a dy- 
namical assignment of the eigenstates wherein some or 
all of the quantum numbers arise by focusing on impor- 
tant quantum and or classical dynamical structures. A 
consistent picture that is emerging from many studies 
is that the complicated spectral splittings and patterns 
at high energies can be 'unzipped' to some extent using 
the underlying classical dynamics. More specifically, the 
spectra are unzipped by recognizing eigenstates showing 
similar localization characteristics about important clas- 
sical invariant structures in the underlying phase space. 

Studies based on classical-quantum correspondence 
have been successfull y appli ed t o sy stems with two cou- 
pled modes[i[l3,[ll|lllll[ll[lI[I| but are yet to be 
extended to systems with three or more coupled modes 



2 



i.e., multimode systems [T7|. The technical and concep- 
tual difficulties associated with a straightforward general- 
ization of the two mode techniques to multimode systems 
are well understood TE'l. Circumventing the technical dif- 
ficulty, in our opinion, requires utilizing quantum objects 
which are sensitive to the underlying classical mechanics 
but do not rely on visualizing the phase space and eigen- 
states. Considerable work has been done in this direc- 
tion and a common theme underlying such approaches is 
the analysis of the eigenvalues and their variation with 
Hamiltonian parametersO IT9II20II2II2II23I . For in- 
stance the nature of an eigenstate |a) of iJ has been 
studied using the methods of diabatic correlations [T9I l20| 
revealing the existence of formal quantum numbers. The 
basic quantity in these studies is the variation of the 
eigenvalue Ea with a specific coupling strength i.e., 
dEa/dTk- Remarkable correlation of the level variations 
to the phase space nature of the eigenstates had been 
noted by Weissman and Jortner in the context of the 
Henon-Heiles system[23|. Support for the correlation was 
provided recently'255 from a semiclassical viewpoint and 
it was suggested that the parametric variations were sen- 
sitive to the various bifurcations occuring in the classical 
phase space. However strong support for the observed 
correlations have existed[2^ in the literature in terms of 
the classical-quantum correspondence of quantum expec- 
tation values i.e., diagonal matrix elements. 

In order to elucidate the connections we note that the 
Hellman-Feynman theorem 



Pk 



(3) 
(4) 



suggests the diagonal matrix element of the perturbation 
as the fundamental object. P^. represents the polyad, 
associated with the resonance Vk, whose constancy is 
destroyed in the presence of other independent reso- 
nant perturbations. Consequently^^ |25 ^ dominance of 
Vj^°'{t; Pk) at a single Pk implies a highly localized \a) 
in the state space with Pk being an approximate quan- 
tum number. Quantum mechanically if |Vj.""| is large 
then it is expected that the perturbation Vk plays a role 
in determining the nature of \a). On the other hand 
a semiclassical analysis of the expectation value pro- 
vides valuable information on the phase space nature 
of \a). This can be seen by considering the quantity 
Pk{E) = J2a Vk°'S{E — Ea) which is the expectation val- 
ues of Vk weighted by the density of states. Such quanti- 
ties lend themselves to an elegant semiclassical interpre- 
tation in terms of the classical closed orbits in the phase 
space. This is hardly surprising given that the genesis of 
the Gutzwiller periodic orbit quantization^^! idea was 
from a semiclassical analysis of the quantum density of 
states p = J2a ^i^^^a)- As the method of semiclassical 
analysis of Pk{E) is well established in the literature and 
our intention in this work is not to semiclassically evalu- 



ate V"" we will highlight the salient features. In general 
Pk has a smooth and oscillating part with the smooth 
part being independent of the nature of the classical dy- 
namics. The oscillating part is sensitive to the nature 
of the dynamics and can be written down in terms of 
a sum over the closed orbits in the phase space. If the 
phase space is chaotic then the closed orbits are the vari- 
ous periodic orbits whereas for a regular phase space the 
closed orbits are the rational tori. In either case it can be 
shown that the oscillating part depends inversely on 
the determinant of the stability matrix M of the closed 
orbit and directly on the quantity 



Vkp 



1 



T 



dtVk{e{t),i{t)) 



(5) 
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representing the average of the resonant term over one 
period Tp of the closed orbit. 

We emphasise the dependence of pk on the quantities 
M and Vkp for two reasons. Firstly, Vkp clearly under- 
scores the important role played by the classical analog 
of Vk- Further, performing a standard canonical trans- 
formation (I, 6) — > (J, ip, K, x) with (J, ijj) being the slow 
angle and action variables specific to the resonant term 
cos(mfc • 6) we obtain 



Vkp = i- 



dtfkiJ,K) cosV' 



(6) 



p ^0 



Now for a closed orbit in the phase space correspond- 
ing to riifc -9^0 the angles x are fast and can be 
averaged resulting in the actions K(i) « K(0). The 
fixed points for the averaged system are then determined 
by (Jpjipp) = (0,0) and correspond to the closed or- 
bit in the full phase space. Within this averaged view- 
point Vkp ~ ±/fc(Jp,K) with the signs coming from 
ipp = 0, ±7r. Thus it is expected that a maximum in 
|V^""| comes from the localization of the eigenstate \a) 
due to closed orbits associated with Vk . Evidently states 
influenced by a particular closed orbit can be classified 
into a group and identified by patterns in the Vj^" "spec- 
trum" . In mixed phase space regimes, generic to molecu- 
lar systems, periodic orbits with varied stabilities can ex- 
ist and influence the dynamics. In the following sections 
we show the expectation values which are scaled to unit 
variance and zero centered. The scaling is performed to 
remove the deoendencep^ |25| of Vj"'" (x p^''^^ on the 
approximate polyad Pk arising from the localization of 
the eigenstate \a) due to Vk. 

Secondly, bifurcations in the phase space are signalled 
by the vanishing^ of the determinant det|(A/ — 1)|. In 
general bifurcations imply birth of new orbits and/or 
death of old orbits. In the molecular context^ such bi- 
furcations manifest themselves as destruction of the old 
modes and creation of the new modes which influence 
the eigenstates. For instance a normal to local transi- 
tion with varying energy arises due to bifurcations in the 
classical phase space. Thus it is natural to expect that 
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the V"" "spectrum" will exhibit the effects of the various 
bifurcations. 

Recent work from our group has demonstrated the 
usefulness of the above method for understanding the 
highly excited states of a model system .25] and the DCO 
radical^^. In this letter two molecules, CHBrClF and 
C2H2, are chosen as further examples to illustrate the 
method. Specifically, the N — 5 polyad of CHBrClF and 
the [16,0]f+ bending polyad of C2H2 are investigated. 
In section II the highly excited states of CHBrClF are 
analyzed and compared to a recent work[l^. We find 
some discrepancy in the dynamical assignments provided 
earlier and by the level velocity method. In section III 
the pure bending states of C2H2 are investigated and it 
is shown that the level velocity approach is capable of 
identifying the new class of local bending and counter- 
rotation states. Comparing to a recent work|23 we show 
that the level velocity approach is successful in identify- 
ing important eigenstate sequences at such high levels of 
excitation. Section IV concludes. 



II. DYNAMICAL NATURE OF CHBRCLF 
EIGENSTATES IN POLYAD TV = 5 

An effective Hamiltonian for CHBrClF was proposed 
by Beil et al. on the basis of their detailed study of 
the rovibrational spectraj28j| . The CH overtone spectrum 
implicated multiple Fermi resonances involving the pure 
CH stretch (s) and the two (a, h) CH bending vibrations. 
The spectroscopic Hamiltonian has the form: 

H = Ho+ ksaaVsaa + ksbbVsbb + kgabVsab + iVaabb (7) 

The zeroth-order anharmonic part is diagonal in the 
number {vs,Va,Vb) basis 

sab sab sab 

■^0 = 5Z ^J^J ^ ^ XjkVjVk (8) 

j j k 

with the i^j and Xjk being the harmonic frequencies and 
anharmonicities of the three modes respectively. The 
perturbations V are off-diagonal in the {vs,Va,Vb) ba- 
sis and represent resonant coupling of the modes. The 
first three perturbation terms in the Hamiltonian are 
Fermi resonances between the CH stretch and the various 
bend modes whereas the last perturbation is a Darling- 
Dennison resonance between the two bend modes. The 
structure of the effective Hamiltonian implies the exis- 
tence of a constant of the motion N ^ Vs + {va + Ufc)/2 
called as the polyad number. For the various param- 
eter values and form of the resonant operators we re- 
fer the reader to the earlier works l2q|. However, we 
note that the bending mode anharmonicities Xaa, Xbb, Xat 
are small and the resonant coupling strengths are rather 
large. For instance, Xaa —6 cm^^ and the stretching 
mode anharmonicity Xss ~ —65 cm~^ whereas ksbb ~ 113 




N N 

FIG. 1: Resonance expectation values as a function of the 
polyad N = Vs + {va + vb)/2 for (a) lowest and (b) highest 
energy states in a polyad for CHBrClF. The circles, squares 
and triangles represent the expectation values of Vaaa, Vsbt, 
and the Vaabb respectively. The IPRs for the two sets of states 
are shown in (c) and (d). Open circles denote IPRs in the 
zeroth order basis and the solid line, dashed line, dash-dotted 
lines represent IPRs in the Ho + Vsaa, Ho + Vsbt, Ho + Vaabb 
basis respectively. 

cm"^. Such large resonant strengths combined with re- 
duced bend anharmonic constants imply that any analy- 
sis solely based on Hq would be inadequate. 

We also report the inverse participation ratios (IPR) of 
the eigenstates in various basis. The IPR La of an eigen- 
state \a) in a basis \b) is given by La = J^b is 
a measure of the extent of delocalization of an eigenstate 
in a specific basis. A high La indicates localization in 
the state space whereas for a completely delocalizcd state 
La = l/Ns with Ns being the total number of states. 

In this work we focus on the Ng = 36 eigenstates 
belonging to the polyad 5 in the energy range of 

[11361, 13947] cm~^. In order to determine the possibil- 
ity of bifurcations with varying A^ in Fig. ^,b we show 
the relevant expectation values as a function of increasing 
A^ for the lowest and highest energy states respectively 
in a given polyad. The expectation values of Vsab are 
relatively small and hence not shown in the figures. The 
IPRs in the zeroth order basis and various intcgrable sin- 
gle resonant basis are shown in Fig. ^,d. For the lowest 
energy state it is clear that the s55- Fermi resonance plays 
an important role. This is evident from the IPR infor- 
mation as well since the sbb single resonance basis values 
are fairly high for the range of A^ shown. On the other 
hand the highest energy states show the importance of 



the saa-Fermi resonance. The sharp rise in the expec- 
tation value corresponding to Vgaa between A'^ = 4 and 
iV = 6 indicates a change in the nature of the dynamics 
which is confirmed by studying the surface of sections. 
In both instances the dominance of a single {V) and the 
corresponding single resonance basis IPR establishes the 
states to be regular and highly localized. 

In Fig.|21the relevant (V) are shown for the eigenstates 
in = 5. Note that a large |V;f"| implies a localized 
state in the corresponding state space resonance zone and 
also in the phase space about a closed orbit. The nature 
of the closed orbit i.e., stable or unstable is inferred from 
the sign of A look at the partial expectation val- 

ues Vi^°'{Pk) will reveal a single dominant contribution 
leading to the assignment (iV, , v) with i/ denoting an 
excitation index for states with the same Pk- We will 
not give a list of assignments for the states in iV = 5 but 
rather focus on certain eigenstate sequences in order to 
compare to a recent workjl^. At first glance it is appar- 
ent that the sbb resonance is influencing the lower end of 
the polyad whereas the saa resonance is important for the 
higher energy states. Note, however, that with increasing 
energy the saa rapidly gains significance with concomi- 
tant decrease in the sbb influence. From Fig. [21 it is seen 
that around 12020 cm~-^ {Vsaa) ~ (Vstb) and the sequence 
splits into two branches. This indicates a change in the 
underlying phase space and thus localization nature of 
the eigenstates on the two branches should be different. 
This is supported by examining the phase space and the 
alternating assignments provided in a recent workp^. 
At the same time towards the middle of the polyad the 
Darling-Dennison resonance is playing a key role in orga- 
nizing the eigenstates. We plot the {V) versus the energy 
eigenvalues Ea since this immediately reveals the states 
that are possibly involved in avoided crossings. In Fig. [21 
certain sequences have been shown and labeled according 
to a classification done recently. In this classification [l^ 
B represent nonresonant states, and C represent states 
influenced by the Darling-Dennison resonance. From the 
expectation values, however, it is not apparent that the 
first few states are nonresonant. Indeed our analysis sug- 
gests some of these states, for instance the first three 
states, are influenced by the sbb Fermi resonance. That 
this is indeed the case has been confirmed by the method 
of local frequency analysis which shows extensive sbb 
locking. In addition the corresponding IPRs in the sin- 
gle sbb resonance basis are very high suggesting strong 
infiuence by the resonance. Thus state number one is 
assigned as {N = 5, Psbb = 2ws + Vb — 10, Va — 0) rather 
than the earlier assignment of {N = 5,Va — 0, Vb = 10). 

The earlier assignment |0 of Darling-Dennison states 
essentially agrees with our results. Certain strongly lo- 
calized sequence of C states is immediately revealed by 
the expectation values as evident from Fig. [21 It is 
expected that such states are localized about the sta- 
ble periodic orbit associated with the Vaabb which has 
been confirmed by computing the Husimi distribuion 
functions. Our assignment for the state denoted C50 is 
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FIG. 2: Resonance expectation values for the states belonging 
to = 5 for CHBrClF. The horizontal zero-line is shown 
for ease of visualization. The symbols used are identical to 
the ones used in Fig. The approximate energy at which 
large scale chaos sets in the phase space is indicated. Some 
of the eigenstate sequences are shown and labeled as in an 
earlier work[T3l. Note the particularly large values for certain 
states indicating localization in state space and phase space. 
Also note that around 12020 cm~^ the expectation values 
(Vaaa) ~ (Vghb) and the splitting of the 'B' sequence into two 
shown by dash-dotted line. Further details are provided in 
the text. 



{N = 5,Paabb = Va + Vb = 10, Vs = 0) and clearly the 
approximate polyad Paabb can be identified with the lon- 
gitudinal quantum number ni introduced in the earlier 
workjl^. Indeed the maximum magnitude of the ex- 
pectation V^""f, can be estimated as P^abb/^ ^^'^ state 
C50 shows close agreement with this classical estimate. 
The so called transverse quantum number rit is associ- 
ated with the degree of excitation v for a given Paabb- 
In the nt = sequence it is clear from the figure that 
with increasing energy other resonances are coming into 
effect. The state labeled C20, in particular, is strongly 
infiuenced by both the saa and the sbb Fermi resonances 
and hence it is inappropriate to classify them as Darling- 
Dennison states. The sequence corresponding to ti; = 1 
is strongly perturbed by the saa Fermi resonance. 

Finally state number 16 and 30 have been classified in 
the previous work as "chaotic" states (class D) exhibiting 
a mixture of class B and class C states. The present anal- 
ysis supports the above classification for state 16, which 
incidentally has the lowest IPR among all the states, but 
clearly Fig. [21 suggests state 30 to be different. The large 
expectation values of saa with smaller sbb expectation 
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values indicates this state to be influenced by the saa and 
sbb resonances. This state can be nominally assigned as 
{N = 5,PsaQ = 2vs + Va = S.viy = 2). Similarly state 
number 33 was classified as E{D) suggesting dominant 
saa character. However state space, expectation values 
(cf. Fig. El and the various IPRs (integrable single shh 
resonance basis IPR is about 0.7) point to a state with 
sbb character. 



III. NATURE OF THE HIGHLY EXCITED 
BEND STATES IN C2H2 

The pure bending dynamics of acetylene at about 
10000 cm^^ above ground state is known to be quite 
complicated due to the strong interaction between the 
trans and the cis bending normal modes However 
recent work by Jung, Taylor and Jacobson revealed that 
the highly excited bending eigenstates were assignable 
despite strong chaos in the svstempoj. The key to their 
dynamical assignments was the fact that a few periodic 
orbits were organizing the dynamics at such high energies 
and hence influencing the structure of the eigenstates. 
Many sequences of eigenstates exhibiting similar local- 
ization patterns about the periodic orbits were identi- 
fied leading to the unzipping of the spectrum. The pur- 
pose of this section is to show that parametric variations 
can immediately reveal the existence of such progressions 
without the need for the visualization of the phase space 
and/or the eigenstates. 

The effective Hamiltonian appropriate for the study 
of the bend only dynamics up to 15000 cm^^ with an 
accuracy of 1.5 cm"! has the formi^ H = + ^'^"'^ + 
i?™*. The harmonic part = u^V/^ + uj^ij^ and the 
anharmonic part is given by 

fjanh ^ ^ (XjkVjVk+gjkljlk)+ ^ VjklVjVkVl (9) 
j,fc=4,5 jkl=i,5 

with Vi and W5 representing the number of quanta in 
the trans and cis normal modes respectively. The de- 
generate bends further require the vibrational angular 
momenta I4 and ^5. Thus the zeroth order states, eigen- 
states of H^"' + m"^, are denoted by Iv^^vi). How- 
ever if™* = Vddi + Vddii + Vii couples the zeroth or- 
der states via the off-diagonal anharmonic resonances. 
Vddi is a Darling-Dennison resonance leading to ex- 
change of quanta between the two modes at constant I4 
and Z5 whereas Vu is a vibrational /-resonance which leads 
to exchange of vibrational angular momentum only be- 
tween the two modes. Vddii results in both exchange 
of quanta and vibrational angular momentum between 
the two modes. For the fitted parameter values and form 
of the perturbatio n op erators we refer the reader to the 
previous works p9ll30j| . Due to the nature of the resonant 
couplings there are two conserved quantum numbers or 
polyads denoted by Nf, = V4 + V5 and / = ^4 -I- ^5. In ad- 
dition eigenstates respect certain symmetries and hence 




FIG. 3: Expectation values as a function of the bend polyad 
Nb = V4 + V5,l = for (a) lowest and (b) highest energy states 
in a polyad for C2H2. Circles, squares and triangles denote the 
expectation values (Vddi), {Vddii), and (Vu) respectively. 
Note the sharp change around A^;, — 12 and Nb — 14 for the 
lowest and highest states respectively. The IPRs in various 
basis for the two states are shown in (c) and (d) as a function 
of Nb- IPR in the zeroth order basis are denoted by open 
circles and in the various single resonance basis by lines {Hq + 
Vddi by solid line, Ho + Vddii by dashed line, and Ho + Vu 
by dot-dashed hue). 



labelled by ± (parity under say I4 — > — ^4) and g/u (v^ 
even/odd). There are a total of 81 states for the polyad 
Ni, = 16 spanning an energy range of [10239,11255] 
cm~^. Although it is possible to analyze all of the eigen- 
states, in this letter we focus on the subset of eigenstates 
with Nf) = 16 and I = and symmetry class g+ i.e., 
states belonging to [16,0]^"*". Note that there are a to- 
tal of 25 states in [16,0]^+. As in the previous section 
we will be interested in the expectation values (three of 
them for acetylene) V^°". 

To begin with we consider, in analogy with the previ- 
ous section, the possibility of various bifurcations occur- 
ing in C2H2 with varying polyad N},. Earlier studies [TH 
have revealed that the lowest (trans bend) and highest 
(cis bend) energy eigenstates in a given polyad undergo 
a sharp change in character with increasing iVf,. More 
precisely, around iVf, = 12 the lowest state is not the 
usual trans bend and becomes a local bend whilst around 
Nb — 14 the highest state changes character from a cis 
bend to the so called counter-rotator mode. Such changes 
are characterisitc of bifurcations which lead to new types 
of modes (dynamics) and we expect the expectation val- 
ues to be sensitive indicators of these changes. In order to 
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FIG. 4: Expectation values for the states belonging to the 
[16, 0]s+ polyad of C2H2. Symbols used are identical to those 
in Fig. El Open Diamonds mark a sequence of eigenstates rep- 
resenting the family I and open circles indicate the sequence of 
eigenstates representing the family II as in an earlier workl^^. 
The lines are drawn as a guide to the eye. Note some of the 
family I states are |u)"'" states shown in this plot for clarity. 
The approximate energy region over which considerable chaos 
exists in the classical dynamics is also indicated. See text for 
discussions. 



show this in Fig.|3^,b we show the {V) for the lowest and 
the highest energy states as a function of Nt,. As expected 
the {V) very clearly indicate the sharp nature of the bi- 
furcations in agreement with the previous observations. 
In Fig. |3t,d we have also shown the IPRs of the states 
in various basis for comparison. Note that the IPRs also 
indicate the transition. However, the various IPRs in 
the post transition regime take on very small values and 
by definition this indicates that many zeroth order nor- 
mal modes are contributing to the eigenstates and thus 
suggesting highly mixed states. Contrast this with the 
fact that all of the expectations are very large suggesting 
highly localized states. The resolution to these appar- 
ently conflicting observations lies in realizing that new 
types of dynamics i.e., periodic orbits created due to the 
bifurcations are influencing the eigenstates and the ze- 
roth order normal modes are a poor basis to understand 
the new modes. It is important to note that IPR in any 
single resonance (dressed) basis is also an insufficient in- 
dicator of the new modes. The expectation values on the 
other hand are quite sensitive and the nature of the new 
modes can be deciphered with further analysis [l^. 

We now turn to an analysis of the states in [16, 0]^+. In 
Fig.^the expectation values are shown for all the states. 



Based on previous observations it is easy to identify the 
local bend and the counter rotating states at the ener- 
getic minimum and maximum of the polyad. The appear- 
ance of rather complicated behaviour in the middle of the 
polyad is related to the occurence of bifurcations with 
varying energy p9j|. Further understanding can be g ained 
by analyzing the reduced classical Hamiltonian|29l| : 

H{J, ;K) = ifo(J;K) + 2s45/a(J;K)cos7A<, 

+ 2r45/6(J;K)cosV'b (10) 
+ 2<45[/_(J; K) cos(V'a - i^b) + /+(J; K) cos(V'a + V-fc)] 

where J = (J^, Jb) = ((w4 - W5)/4, (/4 - k)/^), K = 
{Ka,Kb) = {{Nb + 2)/4, V4), and % = (rs + 2545)/4. 
The angles ipa and V'b correspond to the DDI and the 
vibrational-? resonances respectively. The angle combi- 
nations tpa i ipb correspond to the DDII resonance. The 
fixed points of the above Hamiltonian i.e., (J, ■0) — (0, 0) 
represent the periodic orbits of the full system. It is 
easy to see that in the angle space the fixed points 
are ijja,b = 0, ±7r. From our discussions we expect the 
functions fa, ft, and {f^,f^) at the fixed points to in- 
dicate the 'dynamical' nature of a specific eigenstate. 
In particular, the signs of the various expectations pro- 
vide information on the nature of eigenstate localiza- 
tion in the (ipaTipb) space. It is easy to see that four 

possible sign combinations -|--|--|-,-| , h, and 

— I — are possible corresponding to the localization about 
(0, 0), (0, tt), (tt, tt), and (tt, 0) respectively. At this stage 
we emphasize that highly localized eigenstates imply a 
specific sign combination and large magnitude expecta- 
tion(s). Thus, for instance, from Fig. ^ we anticipate 
that the states \g)f, [5)4 , \g)25, and |g)^y are highly lo- 
calized states about (0, 0), (0, tt), (tt, tt), and (tt, 0) respec- 
tively. These states exhibit localization in phase space as 
well. We find a total of 18 states that belong to one of 
the four sign combinations. On the other hand there are 
states that do not belong to one of the above four classes. 
For example \g)i4^ has the sign combination — h + and 
small expectation values. This suggests a fairly delocal- 
ized state which is hard to assign dynamically. Similarly 

state [5)^9 comes with the sign combination but 

Fig. 21 clearly implies a localized state influenced by the 
vibrational-? resonance alone. Interestingly such states 
are precisely the ones that were multiply assigned in the 
previous studvp^. 

Jung et. al. have used the semiclassical representation 
of the eigenstates [sij in the (^g, il^t) space to provide a 
sequence of similarly localized states |29|. Two main fam- 
ilies, I and II, were identified based on the excitation 
along specific periodic orbits. In Fig.^lwe show two such 
sequences, one from each family, as seen from the per- 
spective of the expectation values. It is apparent that 
the family I is more robust than the family II sequence. 
This is directly related to the fact that the periodic orbit 
underlying family I undergoes far fewer bifurcations as 
compared to the periodic orbit corresponding to family 
II. However even in the family I sequence two states |w)^2 



7 



and |.g)^3 are significantly perturbed. There can be many 
sources for sucfi perturbations and in this case it happens 
to be an avoided crossing. Nevertheless from Fig. 01 it is 
clear that the state \g)fj is localized mainly due to the 
DDII resonance. 



IV. CONCLUSIONS 

In this letter we have demonstrated the utility of eigen- 
state expectation values of resonant perturbations in 
understanding the nature of highly excited vibrational 
states. A particularly large expectation value implies 
that the associated eigenstate is localized and influenced 
by a closed orbit in the classical phase space. Existence of 
localized states and their sequences can be ascertained by 
inspecting the expectation values (magnitude and sign) 
and the IPRs irrespective of wether the classical phase 
space is (near)-integrable, mixed or chaotic. In case of 
complete chaos, from random matrix theory arguments, 
a typical expectation value is expected^^l to be l/^/Wl. 
The examples studied here are quite far from such a 
limit. Combined with the relative signs of the expec- 
tation values it is possible to identify specific eigenstate 
sequences. Sensitivity of the expectation values to the bi- 
furcations provides information on the birth of new types 
of modes and the resulting perturbations on the eigen- 
state sequences. Currently it is not possible to identify 
the type of bifurcation that occurs and this aspect needs 
further study. It is important to note that we are not 
advocating the use of periodic orbits to compute the ex- 
pectation values which is a difficult task. Instead we are 
emphasizing the manifestation and utility of such classi- 
cal structures in a fundamental quantum object i.e., an 



expectation value. Similar philosophy has been adopted 
in the 'vibrogram' or {E, t) approaches to study reso- 
nant svstems|3^ This classical-quantum correspon- 
dence aspect of the expectation values can be easily ap- 
plied to multidimensional coupled systems without the 
need for determining/visualization of phase space, peri- 
odic orbits or eigenstates. Indeed preliminary work on a 
coupled 4-mode effective Hamiltonian suggests that such 
an approach is useful. 

Finally we note that Jacobson and Field have 
recently 35] studied expectation values of resonance op- 
erators for a nonstationary state. Choosing the non- 
stationary state to be a zeroth order bright state |v) 
it was shown that the time-dependent expectations 
(v(t)|V/c|v(i)) indicate the resonances important for dy- 
namics over a particular time interval. Relation to the 
present work is realized by the fact that 

lim i f{^{t)\%Ht)) = ^ |(v|a)pyr (11) 

Thus the time-dependent expectations are, in the infinite 
time limit, nothing but intensity weighted eigenstate ex- 
pectation values. Such an object, the intensity-velocity 
correlator, has been recently introduced and studied in 
great detail|36|. 
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